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The neutrino ground state in a neutron star 
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We address a recent claim that the stability of neutron stars implies a lower bound on the mass of the neutrino. 
We argue that the result obtained by some previous authors is due to an improper summation of an infrared- 
sensitive series and that a non-perturbative "resummation" of the series yields a finite and well-behaved result. 
The stability of neutron stars thus gives no lower bound on the mass of the neutrino. 
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1. INTRODUCTION 

In this talk we present a calculation of the in- 
teraction energy due to multi-body neutrino ex- 
change in a neutron star (][ . 

Consider the series shown in Fig. |l|, which rep- 
resents the "self-energy" of a neutron star due to 
neutrino exchange. In this figure the crosses rep- 
resent insertions of the neutron density. It is not 
difficult to see that the term in this series with k 
insertions of the neutron density scales approxi- 
mately as [|| 
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where N ~ 10 is the total number of neutrons in 
the star, R ~ 10 km is the radius of the star and 
Ck is a dimensionless numerical coefficient. The 
thing which is perhaps surprising in this expan- 
sion is that the "expansion parameter" GfN/R 2 
is of order 10 12 . A direct summation of the terms 
in this series (but truncated after N terms) yields 
an enormous value for the interaction energy . 
This result has led some previous authors to claim 
that neutrinos must have a mass of at least 0.4 eV 
in order to allow neutron stars to exist as stable 
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Our approach to this problem is quite differ- 
ent. We contend that the series represented in 
Fig. |l] is actually infrared-divergent and must be 
"resummed" non-perturbatively in order to yield 
a sensible result. This approach was previously 
advocated in Ref. || . We have performed the re- 
quired resummation and find that the apparent 
infrared divergence is an artifact of the expan- 
sion and that the interaction energy is finite and 
well-behaved |J. There is thus no lower bound 
on the mass of the neutrino. Along the way 
we also encounter some interesting physics. We 
demonstrate, for example, that the ground state 
of the system actually contains a non-zero neu- 
trino number - a result which was previously an- 
ticipated in Ref. 0. In our simple model for the 
density of the neutron star it is straightforward to 
calculate both the energy and the neutrino num- 
ber of the ground state. 

Before describing the calculation in more de- 
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Figure 1. Perturbative expansion of the shift 
in the neutrino ground state energy due to the 
presence of a neutron star. Solid lines represent 
neutrino propagators and crosses represent inser- 
tions of the neutron density. 
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tail, let us note that there have been several 
groups which have examined various aspects of 
this problem |||||--|l2||. In particular, Arafune and 
Mimura |^2| have confirmed our asymptotic result 
using an analytical approximation. 

2. NEUTRINO GROUND STATE 

2.1. Preliminaries 

Our goal is to calculate the shift in the neutrino 
ground state energy due to the presence of the 
star. This energy shift may be defined in terms 
of the neutrino Hamiltonian H( ) in the presence 
(absence) of the star as follows |13| : 

W = (6\H\6) - (0\H \0). (2) 

Here |0) denotes the neutrino ground state in the 
presence of the star, while |0) denotes the usual 
matter-free vacuum state. As we have already al- 
luded, the state |0) contains in general a non-zero 
neutrino number (i.e., it is "charged"). Note that 
the expression in Eq. (|J) is a formal, ultraviolet- 
divergent quantity which needs to be renormal- 
ized. This renormalization may be done quite 
easily using the usual techniques. 

In order to proceed, it is convenient to in- 
troduce an effective Lagrangian for the neutrino 
field. After integrating out all of the other parti- 
cles in the theory, one obtains 

C eS =^ L [ty+arf] Tp L (3) 

where ipL — ^(1 — Ts)^) an d where 

a(x) = G FPn {x)/V2 ~ 20 eV (4) 

is the electroweak potential induced by the fi- 
nite neutron density (p n ss 0.4 fm~ 3 in a typ- 
ical neutron star). This potential is iden- 
tical to the one which is usually considered 
in the well-known Mikheyev-Smirnov-Wolfenstein 
(MSW) effect |§. Note that the potential term 
in (|^) resembles a position-dependent chemical 
potential, so that it is not at all surprising that 
the ground state of the system has a non-zero 
neutrino number. 

It is straightforward to derive the follow- 
ing perturbative expansion for W in terms of 
the potential a(x) and the neutrino propagator 
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W = -g-y^^IoGbrf (5) 

= J2 wik} - ( 6 ) 

fc 

All of the odd terms in this expansion disappear 
so that this series corresponds precisely to that 
which is represented diagrammatically in Fig. ^. 
One aspect of the perturbative expansion which 
is useful is that it neatly isolates the ultraviolet 
divergence in W . In fact, the only ultraviolet di- 
vergent term in (||) is that with k — 2. This 
term is related to the vacuum polarization of the 
Z in the complete theory. A final note concern- 
ing this expansion is that while the terms with 
k > 4 are separately finite, their sum is infrared- 
sensitive and is in fact ill-defined for "large" stars 
[aR 3> 0(1)]- The non-perturbative "resumma- 
tion" of these terms is the main goal of our cal- 
culation. 

2.2. Comparison point 

It is useful to compare the results which we will 
describe with results which were obtained pre- 
viously in the literature. The perturbative ex- 
pansion of W has been considered in Refs. (| ||]. 
where it was found that already by the eighth 
term in the expansion the interaction energy ex- 
ceeded the gravitational binding energy of the 
star. After summing up the entire series [which in 
the approach of Refs. was actually a trun- 

cated sum, not an infinite series as we have in 
Eq. (^)], it was found that the interaction energy 
exceeded the rest mass of the universe. It was ar- 
gued that the only way to regulate the sum was 
to give all neutrino flavours a mass of at least 0.4 
eV. 

Our philosophy in this matter is that the per- 
turbative expansion is simply outside of its ra- 
dius of convergence when aR ^> 0(1) and that 
the series needs to be resummed using a non- 
perturbative approach ||. We cannot consider a 
realistic neutron star (i.e., with aR ~ GfN/R 2 ~ 
10 12 ) using our numerical approach, but it is ac- 
tually sufficient to restrict our analysis to aR < 
0(100), since for aR ~ 0(1) we already observe 
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a "cross-over" to the non-perturbative regime. A 
clear signal of this cross-over is that the ground 
state obtains a non-zero charge. (The charge of 
the ground state is exactly zero to any finite order 
in perturbation theory.) 

Consider the comparison point aR = 20. This 
point could correspond to a tiny "star" with 
a realistic neutron density (a — 20 eV), but 
with a tiny radius (R = 2 X 10~ 5 cm). In 
this case the truncated sum in Ref. || gives 

J2k=4 w(k) ~ 1066 eV - B y W of compar- 
ison, our non-perturbative resummation gives 

Ek=4 w(k) ~ -2- 3 keV - 

2.3. Phase shift formulas 

An exact, non-perturbative expression for the 
interaction energy W is given by the following 
expression due to Schwinger (l3| 

W = — V (21 + 2) / (kj [Si (w) + Si (-oj)} . (7) 
27T J=t Jo 

Here 8i(u>) is the scattering phase shift for a neu- 
trino incident on the star and I labels the orbital 
angular momentum. A similar expression may be 
obtained for the "charge" of the ground state: 



g = -i ? 2J(2i + 2) 



1=0 



duj 



(8) 



The factor (21 + 2) = (2j + 1) is the degeneracy 
factor for a given energy uj and total angular mo- 
mentum j. 

The beauty of the above expressions for W and 
q are that they are valid for any value of aR. 
Note, however, that Eqs. (0) and (||) are still for- 
mal ultraviolet-divergent expressions which need 
to be renormalized. 

3. NUMERICAL EVALUATION 

Let us first choose a density profile for the neu- 
tron star. Since the quantities which are of in- 
terest to us (i.e., W and q) are only sensitive to 
the gross features of the star, it is convenient to 
choose a very simple - albeit unrealistic - density 



profile: 

a(r) = a6(R-r). (9) 

Recall that a = G F p n (x)/V2 ~ 20 eV in a neu- 
tron star. In this very simple model it is straight- 
forward to obtain closed expressions for the scat- 
tering phase shifts, which simplifies our numerical 
work considerably. 

We may now renormalize W and q. Since 
our model is renormalizable the ultraviolet di- 
vergences in W and q are confined to the first 
few terms in the perturbative expansion. These 
terms may be isolated by Taylor-expanding the 
phase shift formulas in aR. [Note that by Taylor- 
expanding in aR we recover the perturbative ex- 
pansion defined in Eq. (g).] The procedure is then 
as follows: (i) Taylor expand W and q in order to 
isolate the divergent terms; (ii) subtract out the 
divergent terms; (iii) regularize and renormalize 
the divergent terms using conventional methods; 
(iv) add the finite, renormalized terms back in. 
This procedure yields the following expression for 
the renormalized energy: 



Wren = 

where 



w& + W {A+) 

" rcn ' 



(10) 
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— J2(2l + 2) dw[5,(w)+5,(-«)],(ll) 



and where Si(lo) = Si(lu) — 5j 2 \u>), with 8Y'(lj) 
being the second term in the Taylor expansion of 

(2) 

Si(u>). Wren is essentially the vacuum polariza- 
tion of the Z in the full theory (but convoluted 
over the neutron star) and corresponds to the first 
term in the diagrammatic expansion in Fig. 1 . 
This term has been discussed in detail in Ref. 1] 
and will not be considered further here. The sec- 
ond term in Eq. ©, W^ A+ \ is of particular inter- 
est to us since it represents the non-perturbative 
"resummation" of the terms with four or more 
insertions of the neutron density in the diagram- 
matic expansion shown in Fig. [I]. As we have 
noted, a direct summation of these terms does 
not lead to a sensible result when aR > 0(1). 
The resummed result in Eq. (pi]), however, is al- 
ways well-defined and leads to a well-behaved and 
sensible result. 



?(2) 
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We may similarly obtain an expression for the 
renormalized charge: 

oo 
1=0 

In order to calculate the charge, then, we need 
only know the scattering phase shifts at the ori- 
gin. 

3.1. Small aR 

For small aR the "resummed" energy, W (i+ \ 
is well-approximated by the leading term in the 
Born expansion. We have checked this explic- 
itly by calculating W^ 4+ ^ directly using the phase 
shift formula and by Taylor-expanding the phase 
shifts to fourth order and (analytically) integrat- 
ing the resulting expressions. The result which 
we obtain is 

W l<+) ~ W W ~ _o.00115^^. (13) 

R 

For small aR the phase shifts are zero at the ori- 
gin, so that the ground state of the system is un- 
charged [see Eq. (Q)]. 

3.2. Larger aR 

As aR is increased, there is a critical point be- 
yond which it becomes energetically favourable 
for the neutrinos to "condense"; thus, for 
aR> aR\ cHt , the ground state of the system car- 
ries a non-zero neutrino number. There are in 
fact an infinite number of points at which the 
charge of the ground state changes discontinu- 
ously as aR is increased. In our simple model 
for the density profile of the neutron star there 
is a correspondingly simple condition for a new 
charge to be added to the ground state: 

3i(aR) = 0. (14) 

In the quantum mechanical scattering problem, 
these values of aR correspond to the points at 
which a resonance crosses from the positive to 
the negative energy continuum 

Figure || shows a plot of the renormalized 
charge as a function of aR. The solid curve gives 
the exact charge, which has periodic jumps, and 
the dashed curve gives the charge expected in the 
large volume limit for a system with chemical po- 
tential [i = a: <7cond = 2(ai?) 3 /(97r). Clearly, as 
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Figure 2. Plot of the charge as a function of aR. 
The solid curve gives the exact charge and the 
dashed curve gives the charge expected for a sys- 
tem with chemical potential fi — a. 



aR gets large the exact result tends to this limit. 
The critical point - at which the ground state be- 
comes charged and beyond which we should not 
trust the perturbative expansion - is seen to oc- 
cur at aR = 7r in this model. This is the first zero 
of j (aR). 

Figure || shows plots of both the charge and the 
energy as functions of aR. Both the charge and 
the energy are normalized to the values expected 
for a "condensate" in the large- volume limit. For 
small aR the charge is zero and the energy is 
well-described by the leading term in the Born 
expansion [which is shown by the lower dashed 
line in Fig. |^(b)]. At aR — tt the ground state 
becomes charged and perturbation theory breaks 
down. As aR is increased, both the charge and 
the energy tend toward the values expected for 
a large system with chemical potential fi = a. 
[Note that W cond = -a 4 R 3 /(I8ir).} The asymp- 
totic trend apparent in Fig. |3j has recently been 
confirmed by Arafune and Mimura using an ana- 
lytical approximation [ fL2"f . 

The comparison point which was singled out 
above (aR = 20) is included in the points shown 
in Fig. H A previous (perturbative) calculation 
of the energy associated with this point yielded 
an enormous value [9 , while our non-perturbative 
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Figure 3. Normalized plots of (a) the charge and 
(b) the energy as a function of aR. The dots 
give the results of our exact (non-perturbative) 
calculations. 



approach yields a very small and innocuous value. 

4. CONCLUSIONS 

It has been argued in Refs. that the in- 

teraction energy due to the exchange of massless 
or very light neutrinos in a neutron star would 
be enormous and would destabilize the star. We 
have addressed this claim by performing an ex- 
plicit non-perturbative calculation of the inter- 
action energy. We find that, once properly re- 
summed, the interaction energy is a finite and 
well-behaved quantity which is far too small to 
have any effect on the fate of a neutron star. 



Thus, contrary to previous claims, the stability 
of neutron stars places no lower bound on the 
mass of the neutrino. 
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